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ABSTRACT 

Using a large sample of galaxies from the SDSS-DR7, we have analysed the alignment of disk galax- 
ies around cosmic voids. We have constructed a complete sample of cosmic voids (devoid of galaxies 
brighter than Mr - 5 \ogh - -20.17) with radii larger than 10 Mpc up to redshift 0.12. Disk galaxies 
in shells around these voids have been used to look for particular alignments between the angular momen- 
tum of the galaxies and the radial direction of the voids. We find that disk galaxies around voids larger 
than > 15/1"' Mpc within distances not much larger than 5 Mpc from the surface of the voids present 
a significant tendency to have their angular momenta aligned with the void's radial direction with a sig- 
nificance > 98.8% against the null hypothesis. The strenght of this alignment is dependent on the void's 
radius and for voids with < \5h^^ Mpc the distribution of the orientation of the galaxies is compatible 
with a random distribution. Finally, we find that this trend observed in the alignment of galaxies is similar 
to the one observed for the minor axis of dark matter halos around cosmic voids found in cosmological 
simulations, suggesting a possible link in the evolution of both components. 

Subject headings: Large Scale Structure: Voids; Galaxies : General 



1. Introduction 

The study of the alignment of galaxies with respect 
to the large scale structure is a recurrent topic still not 
fully settled. The first works studying the alignment of 
galaxies focused on clusters and superclusters. Stud- 
ies on the alignment of galaxies in clusters (Adams 



et al.,1 980) and superclusters (Flin & Godlowski 1986| 
Kashikawa & Okamura|p^92 1 claimed to find partic 



ns] 

;6f 



ular alignments of the galaxies with respect to their 
local large scale structure. On the other side, similar 



studies did not find any particular aUgnment (Helou & 



|Salpeter|1982[|Dekel|1985t|Garrido et al.|1993| l. More 



recently, Navarro et al.' (2004') revisited the analysis 
done by Flin & Godlowski (1986) on the alignment 
of galaxies in the Local Supercluster (LSC) under the 
light of the Tidal Torque Theory (for a recent review 
about the Tidal Torque Theory or TTT, see Schafer] 
20091. The authors found a tendency of galaxies to 



have their spin parallel to the plane of the LSC, also 
known as supergalactic plane, that would support the 
predictions from the TTT. 

However, the observational analysis is hindered by 
two main difficulties: the accurate determination of the 
direction of the angular momentum of the galaxies and 
the determination of the distribution of matter around 
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them. The determination of the spin of disk galaxies 
can be guessed by the shape of the galaxy, considering 
that galaxies spin around their minor axis. However, 
there is still an indetermination due to projection ef- 
fects since in most of the cases it is not possible to 
know which half, of the two in which a galaxy is di- 
vided by its major axis, is closer to the observer The 
presence of dust lanes or the use of kinematic data can 
help to solve this degeneracy but in most of the cases 
this information is not available. To deal with this 
problem some authors have taken all the possibilities 
of the spin as independent ones ( [Kashikawa & Oka^ 



1992 1 while others have opted for taking just one 



possibility ( Lee & Erdogdu|2007) l. 

Regarding the accurate determination of the mass 
distribution around the galaxies, the main problem 
comes from the effects of the proper motion of galaxies 
which introduces uncertainties in the conversion from 
redshift to distances. 



To overcome both problems, Trujillo et al. (2006 



hereafter T06) proposed the use of spiral galaxies seen 
edge-on or face-on (so the direction of the spin vec- 
tor is better determined) located in the shells around 
cosmic voids. The advantage of the regions around 
large cosmic voids is that the direction of the gradient 
of density is strongly aligned with the radial direction 
which can be determined in a robust way despite the 
uncertainties of converting redshifts in distances. Us- 
ing this technique and data from the third data release 
of the Sloan Digital Sky Survey (SDSS-DR3) and the 



2dF Galaxy Redshift Survey (2dFGRS), |T06| found a 
tendency of galaxies around shells of voids to have 
their spin vector perpendicular to the radial direction. 
Cuesta et al. ( 2008| l working on cosmological simula- 



tions of dark matter halos around voids found results 
in apparent agreement with those of T06, The simu- 
lations show that the angular momentum of the dark 
matter halos tend to be also aligned to the perpendicu- 
lar direction. In both cases, the results were in agree- 



ment with the prediction done using the TTT ( Lee & 



|Pen"2000), that the angular momentum would tend to 
be aligned with the intermediate axis of the tidal shear 
tensor, that in the surface of the voids is in the perpen- 
dicular direction. 

However, recently, |Slosar & White ( 2009 hereafter 
S09) have redone a similar analysis, but using a larger 
sample of galaxies from the SDSS-DR6, obtaining a 
result that is consistent with a random distribution of 
orientations, in contrast with the previous results. 

In this work, we revisit the analysis of the aUgnment 



of galaxies around voids with two significant improve- 
ments with respect to those two previous works. First, 
we make use of the latest data release of the SDSS, i.e. 
SDSS-DR7, and we combine it with the morphologi- 



cal classification from the Galaxy Zoo project (Banerji 



et al.|2010 Lintott et al.[2010|l to select disk galaxies. 



Second, we have developed a statistical procedure to 
partially correct the indetermination in the spin direc- 
tion due to the projection eff'ect so we can obtain infor- 
mation also from galaxies that are not edge-on or face- 
on, increasing by a factor of 3 the effective number of 
galaxies that are used in our analysis with respect to 
the restriction to edge-on and face-on galaxies. 

The outline of this Paper is as follows. Section |2] 
presents the data used for our analysis; Section l3] de- 
scribes the procedure to search for voids; Sectionfflis 
devoted to the selection of galaxies and the computa- 
tion of their alignments; Section [5] contains the final 
results; in Section|6]the results are discussed and com- 
pared with previous works and in Section |7] the sum- 
mary of the resuls are presented. 

Through this paper we assume a ACDM cos- 
mological model with Qm - 0.3, Qa = 0.7 and 
Ho = lOO/ikms-'Mpc-i. 

2. The data 

On what follows we describe the data that we have 
used to: a) create a sample of cosmic voids, and b) 
obtain a sample of galaxies in the shells surrounding 
them to explore the orientation of galaxies in the large 
scale structure. 

Our main source of data has been the New York 
University Value-Added Galaxy Catalo J^NYU-VACG, 



Blanton et al. 2005 ; Padmanabhan et a l.|2008HAdelmarj - 
^McCarthy et al. 2008 ), which is based on the photo- 
metric and spectroscopic catalog of the SDSS-DR'^([sJrauss 
|et al.|p002] l. The main characteristics of the NYU- 
VACG are: 

• Spectroscopically complete up to r ~ 17.8 (ex- 
tinction corrected) 

- Completeness ~ 99% 

- Success rate ~ 99.9% 

• yU5o('' - band) < 24.5mag arcsec"^ 



http://sdss.physics.nyu.edu/vagc/ 
^http://cas. sdss.org/astrodr7/en/ 
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• ~ 90 targets/deg^ 

• Median(z)=0.104 

2.1. Selection of the sample of galaxies 

We have established thresholds in absolute magni- 
tude, M'"", and redshift, z''™, in order to maximize 
the amount of galaxies while keeping the final sam- 
ple complete in spectroscopy. The completeness of the 
initial catalog is a basic requirement to avoid the de- 
tection of spurious voids. 

In Figure[T]it is plotted the number of galaxies with 
Mr < M';'"(z) as a function of the redshift. The value of 
MI.""(z) corresponds to completeness limit (r = 17.8) 
at each redshift z- From the peak of this distribution 
we obtain the limits of our final sample: 

• z<0.12 

• -5 log/i < -20.i:Q 

2.2. Volume trimming 

Although we focus our analysis in the largest con- 
tinuous volume of SDSS-DR7, the irregular limits of 
the volume still posed difficulties in the reliable detec- 
tion of voids. For this reason, we have defined new reg- 
ular limits minimizing the detection of spurious voids 
while still keeping ~ 90% of the original volume. 

Figure |2] shows the limits of our trimmed volume, 
projected onto the original distribution of galaxies, de- 
fined as follows: 

• 6>0° [Southern Umit] 

• S< -2.555556 (a - 131°) [Western limit] 

• 6< 1.70909 (a - 235°) [Eastern limit] 

0.93232 sin(Q-95.9°) 



6 < arcsin 
limit] 



Vl-[0-93232 cos(q--95.9°)]- 



[Northern 



Table [T] summarizes the main properties of the ref- 
erence catalog that we have used in our search for cos- 
mic voids. 



^We opted for no applying a k-correction due to the small redshift 
range probe and the high uncertainties in its determination. 



2.3. Homogeinity check 

Another important requirement of the galaxy cata- 
log to be suitable to search for voids is the homogein- 
ity. One common test of homogeinity is the {V/V„iax} 
test, for which a value of 0.5 is expected for an ho- 
mogenous distribution. 

To perform this test, first of all, for each galaxy it is 
computed the volume, V, of the sphere with radius the 
distance along the line-of-sight to it. Then, the max- 
imum of all the volumes, Vmax, is found and the ratio 
V/V„,ax is obtained for each galaxy. The final step is to 
calculate the average value of these ratios, {V/V,„a^}. 

For our galaxy catalog (V/Vmax) = 0.49990, con- 
firming that at large scale the distribution of galaxies 
in our volume is homogeneous. 

3. Catalog of voids 

With the catalog of galaxies described in the previ- 
ous section, we proceed to search for cosmic voids on 
it. In this section we describe the procedure followed 
to construct our catalog of voids. 

3.1. Procedure description 

First of all, we need to establish the definition of 
"void" that we use in our analysis. We have opted 
for the simplest one: a spherical volume devoid of 
any galaxy brighter than our completeness limit. This 
definition has been already used in other works such 
asJPatiri et al .j (|2006a|l, [Patiri et al.| ( |2006b|l, [Trujillo 
let ai.| ( |2006( ), |Brunino et al.| ( |2007| ), and |Cuesta et al. 
( |2008| l. ICuesta et aL] ( |2008| ) found that for dark matter 
haloes in cosmological simulations, using ellipsoidal 
voids instead of spherical ones does not affect signifi- 
cantly their results. This gives us confidence in the use 
of spherical voids for our analysis. 

Apart from minor differences, the procedure that 
we have followed is basically the HB Void Finder de- 
scribed in [Patiri etaL] ( |2006a| l. These are the basic 
steps: 

1 . Random points are thrown within the volume of 
the catalog. 

2. For each trial point, the 4 closest galaxies are 
found and the center and radius of the sphere 
defined by these 4 galaxies are computed and 
stored. 

3. Of the resulting spheres, those fulfilling any of 
the following criteria are rejected: 
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Reference catalog 


NYU-VAGC (Galaxies) 


Spectroscopic completeness limit 


r< 17.8 


Redshift limits 


0.005 <z< 0.12 


Absolute magnitude limit 


Mr -5 log/i < -20.17 


Number of galaxies 


142127 


Total Projected Area 


1.941484 strad 




0.1545 x47r 


Total Volume 


0.0276556 (/i"' Gpc)^ 


Average density of galaxies 


0.00514 (r'Mpcr^ 



Table 1: Summary of properties of the reference catalog. 



-18.16-19.74-20.69-21.38-21.93-22.38-22.77 




KOO 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
Z 



Fig. 1. — Number of galaxies brighter than M'"" with 
z < z""\ M';"' is the absolute magnitude corresponding 
to the spectroscopic limit (r - 17.8) at redshift z''"'. 
The peak of the distribution it is used to establish the 
thresholds in Mr and z of our initial catalog. 




100 120 140 160 180 200 220 240 260 280 
aC) 



Fig. 2. — Partial projection of galaxies from the SDSS- 
NYU catalog. The Umits used in this work are over- 
plotted. 



• Not being empty. 

• Intersecting the border of the volume]^ 

• Having a radius smaller than 10/; ' Mpc. 

Finally, to have well defined spherical voids, we im- 
pose that they can not overlap. In the case in which 
several voids overlap, only the largest one is kept. The 
order in which the rejection of the voids is done affects 
their final sample. Therefore, to ensure that our final 
sample contains the largest possible voids, the process 
of rejection of overlapping voids it is done from the 
largest void to the smallest one. 



'^Patiri et al.'(2006a^ used a different technique. They put artificial 
galaxies in the limits of the survey, allowing voids to be defined by 
three real galaxies and one fake one. 
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In summary, the voids selected for our final catalogs 
fulfill the following conditions: 

1 . They are empty of galaxies from the initial cata- 
log, i.e. M,. < -20.17 + 5 \ogh. 

2. Their radius is larger than lO/z"' Mpc. 

3. They are completely inside the surveyed vol- 
ume 

4. Voids do not overlap, i.e. the distance between 
the centers of two voids is larger than the sum of 
their radii. 

The power of this procedure to produce a complete 
catalog of voids depends critically on the relation be- 
tween the density of galaxies, the size of the voids 
and the number of trial points used in the search. We 
have performed several tests and have found that using 
~ 10^ initial random points (corresponding to a den- 
sity of ~ 35(/;"' Mpc)"^ trial points) ensures that the 
completeness of our catalog of voids is > 99%. 

The final outcome of this procedure is a complete 
catalog of 699 non-overlapping cosmic voids with ra- 
dius larger than 10 /i"' Mpc. We found a median radius 
of 11.85/1"' Mpc and the average density of voids is 
32.2 X (/j-'Mpcr^j^ 

Table |2] provides the main information of these 
voids. For each void, we include the position of 
the center, both in Cartesian coordinates {X, Y, Z), and 
equatorial coordinates {a, 6) and redshift z; and the ra- 
dius {R). The Cartesian coordinates are computed as 
follows; 



X = 
Y = 
Z = 



D(z) cos 5 cos a 
D{z) cos S sin a 
D{z) sin 5 



where D{z) is the comoving radial distance. 

4. Catalog of galaxies in shells around voids 

Using the previous catalog of voids, we extract 
those galaxies within shells of 10 /i ' Mpc around each 
of them. The morphology of the galaxies has been 



^Be aware that this criterion impUes that the effective volume in 
which the center of voids can reside is smaller than the whole vol- 
ume and depends on the sizes of the voids. 

^Density computed within the effective volume in which the centers 
of the voids can be located. 



obtained from the Galaxy Zoo catalog ( Banerji et al. 
|2010[ [Lintott et al.|201 0l which provides robust dis- 
tinction between elliptical and disk galaxies, although 
~ 52% of the galaxies remain classified as "uncertain". 
We select only those galaxies classified as "spiral" in 
the Galaxy Zoo catalog (~ 37%) because for them the 
direction of the angular momentum is well defined by 
their minor axes. 

This results in a final sample of 45522 measure- 
ments, from 32374 single galaxies, that has been used 
to study the alignment of their angular momentum 
with respect to the void|^ The use of measurements of 
galaxies falling in more than one shell is justified be- 
cause the small uncertainties introduced are compen- 
sated by the increase in the size of the sample. 

In the following sections we describe in detail how 
this alignment has been computed. 

4.1. Computation of the alignment 

The advantage of studying galaxies around voids is 
that, on average, the density increases radially. This 
makes the radial direction a good proxy for the distri- 
bution of matter around each galaxy. Therefore, we 
use the minor axis of the galaxies to define the orien- 
tation of their angular momentum, and the radial di- 
rection of the voids to characterize the distribution of 
matter around them. Hence, our analysis is focused on 
the angle between these two directions, 0^ For prac- 
tical reasons, in the context of this work we will use 
the expression "radial direction" to mean the direction 
defined by a galaxy and the center of the correspond- 
ing void, and "perpendicular direction" to mean any 
direction perpendicular to the radial direction. 

To compute the angle we need to define the di- 
rection of the angular momentum or spin, s, of each 
galaxy. We do this, first of all, by computing the incli- 
nation angle between the plane of the galaxy disk and 



the lin e of sight, Following Haynes & Giovanelli 



( 1984 1 and |Lee & Erdogdu ( 2007| | we use a model of 
thick disk with a projected minor- to-major axis ratio 
a/b and an intrinsic flatness / (i.e. the ratio between 
the real minor axis and the real major axis). According 



^These data is avalible upon request to the authors. 

schematic illustratio n of this method is shown in Figure 1 of T06 



Lee & Erdogdu 



2007 



used the angle i defined as the angle between 
the plane of the galaxy and the projected plane of the sky, therefore 
f = ;r/2 - 
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X 

Mpc] 



Y 

[/!-' Mpc] 



Z 

[r 1 Mpc] 



R 

[h-^ Mpc] 



a 
[deg] 



5 

[deg] 



-193.463 
-227.081 
-226.366 
-86.105 



-145.544 
26.372 
17.375 
93.444 



193.384 
140.238 
175.425 
136.020 



18.703 
18.205 
18.104 
17.606 



216.954492 
173.375602 
175.610699 
132.659407 



38.617468 
31.526741 
37.692842 
46.949188 



0.1059 
0.0914 
0.0979 
0.0630 



Table 2: Catalog of voids (excerpt). Complete version in electronic form. See text for more details. 



to this model, ^ can be obtained applying the formula: 

(b/af-f 



sin^ f = 



(1) 



For values of b/a < f, the angle is set to 0. 

The flatness / depends on the morphological type 
of the galaxies and we use an average value of 0.14. 

From Equation ([T]i it is easy to see that to a single 
value of b/a corresponds two values of = \(\ 

and - -\(\. The indetermination is irrelevant for 
^ - (edge-on galaxies) and for ^ = +n/2 (face-on 
galaxies). 

To compute the spin vector s we followed the pre- 
scription by 'T06' According to this, if {a,d) are the 
equatorial coordinates of a galaxy, ^ the inclination an- 
gle obtained from Equation ^ and <p is the position 
angle of the galaxy increasing counterclockwise (i.e. 
from north to east in the plane of the sky), the compo- 
nents of s are: 



S y- - COS a cos 6 sin f 

H- cos ^(sin cos a sin (5 - cos sin a) (2) 
Sy = sin a cos S sin ( 

+ cos ^(sin (p sin a sin 6 + cos cos a) (3) 
Sj; = sin (5 sin ^ - cos ^ sin (p cos 6 (4) 

Next, we compute the angle between the radial vec- 
tor that connects the center of the void, r^jv/, with the 
center of the galaxy, rgaiaxy'- 



'* ~ f galaxy ^void- 



(5) 



Having obtained r and s, the angle between them, 0, 
it is computed as: 



6 - arccos 



s ■ r 



(6) 



4.2. Analytical model of the distribution of 9 

We compute the angle 9 for all the galaxies in our 
sample of galaxies around voids, obtaining a distribu- 
tion of 9. |Betancort-Rijo & Trujillo| ( |2009| l provide an 
analytical model for the distribution of the angle 9, or 
to be more precise, of |cos 9\, f (|cos 9\). From theoret- 
ical principles confirmed by simulations, the authors 
found that P(|cos 6*1) is well described by the expres- 
sion: 



Pip) 



pdu 



[1 -H(p2- l)^2]3/2 



fi = |cos 9\ , 



(7) 



where p is a free parameter that describes the overall 
shape of the probability distribution. An interesting 
property of this distribution is the relation between the 
parameter p and the average value of |cos 9\: 



1 



<|COS0|> 



1 



(8) 



Besides, the values of p are related with the exis- 
tence or absence of particular alignment between the 
vectors r and s according to the following criteria: 

p < I. r and s tend to be parallel. 

p - I. There is no particular alignment between r and 
s. 

p > I. r and s tend to be perpendicular. 

It has been found that Equation (j7]i describes well 
the results from cosmological simulations ( |Brunino| 
let al.|2007t|Cuesta et al.|2008] l. 

An alternative expression for the same distribution 
described by Equation (j7]i it is provided by Lee ( 2004 1 



and used in several works (eg., T06 Lee & Erdogdu 



2007, |S09| l. This alternative expression is character- 
ized by a parameter c and, from the comparison be- 
tween Equation (9) of S09 and Equation (|7| of the 
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present work, it is possible to obtain the following ex- 
pression relating both characteristic parameters: 



3c 



2(1 - c) 



(9) 



4.3. Statistical computation of P(|cos 611) 



Different approaches have been used to deal with 
the indetermination of the values of ^. For example, 
[Kashikawa & Okamura| ( [T992| i uses the two values of 
^ independently. Another possibility is to use just one 
sign in the definition of ^ as done by |Lee & Erdogdu| 
( 2007[ l. However, these authors acknowledge that this 
decreases the strength of the measured alignment. 

On the other hand, |T06l and |S09| have overcome the 
problem with the indetermination of the values of ^ us- 
ing only edge-on and face-on galaxies, for which the 
direction of the spin is well determined. The main dis- 
advantage of this approach is that the number of galax- 
ies suitable for computing P(|cos 6\) is greatly reduced. 
For example, using the criteria from |T06 the fraction 
of galaxies that can be used is ~ 22% of all the disk 
galaxies. 

We opted for a statistical approach that allows to 
compute a corrected distribution P^dcos d\) from the 
combination of the distributions P(|cos 0\) obtained us- 
ing both signs. 

In Appendix|A]we describe in detail this procedure. 
We also show in this appendix the results from several 
Monte Carlo simulations that show the ability of the 
procedure to recover the correct values of p (Table [7]|. 

The fact that we actually do not know the real val- 
ues of ( is reflected in the uncertainties of the proce- 
dure. Although the values of p are well recovered, the 
uncertainties measured from the simulations are larger 
than those expected from considering just the size of 
the sample, A^^. Of course, this is because we are not 
using the real values of Nevertheless, from the sim- 
ulations we have obtained that our procedure has a pre- 
dictibility power equivalent to that of a sample 0.6 Ng 
with complete knowledge of the real values of f . Let's 
remember that the common procedure of using only 
face-on or edge-on galaxies is restricted to ~ 20% of 
the total amount of spiral galaxies. This means that 
our statistical procedure increases by a factor of 3 the 
effective number of galaxies with respect to previous 
works. 



5. Results 

Using the procedure described in Appendix [A] we 
have computed the corrected distribution Pc{\cos0\) 
of the sample of galaxies around voids. Given the 
large size of our initial sample, we have also computed 
f c(|cos d\) for different subsamples combining differ- 
ent sizes of voids and shells around them. 

An important point of our analysis has been to es- 
tablish the significance of our results in a robust way. 
This is done by comparing, in each case, the mea- 
sured signal with the standard deviation cr^icosei) of 
the theoretical distribution in case of null signal, i.e. 
{|cos0|) = 0.5{^ In the situation of complete knowl- 
edge of the real values of we would have cr(|cose|> = 
( -^12 X Ngy^ for a sample of Ng measures. However, 
we have akeady shown that our procedure has uncer- 
tainties equivalent to a sample of size 0.6 A^^, therefore, 
the previous expression needs to be corrected to 



O" (|cosfl|> 



7l2x0.6xA^g 



(10) 



Knowing the value of cr(|(,ose|>, we can establish the 
signal to noise ratio (SNR) of the signal of a subsample 
of measurements as: 



SNR 



O.5-<|cos0|),, 

O" (|cos fl|> 



(11) 



where {\cos 6\)corr is obtained from the statistical 
correction and cr^icosfli) from Equation ( 10 1. Note that 



the denominator is the signal, which corresponds to the 
difference between the observed value of (|cos 6\)corr 
and that of the random distribution which is 0.5. For 
practical reasons, the sign of SNR has been chosen so 
that is positive for values of /9 > 1 ((|cos ODcrr < 0.5) 
and negative for values of p < I ((|cos 6\}corr > 0.5). 

In Tables |4] and |5] are shown the main results of our 
analysis. In Table |4] samples are constructed by set- 
ting a minimum value for the radii of the voids {R'y"ij) 
while, in Table |5] samples are constructed using voids 
with radii in the ranges Rvoid ± 0.5 /j"' Mpcp] Apart 
from this difference in the definition of the first col- 
umn, both tables share the description of the rest of the 



'"in the analysis of the significance is more convenience the use of 
(|cos f|> than that of p because the fonner has a gaussian distribution 
but the latter has not. 

' ' For convenience, we will refer the samples of the first table as "cu- 
mulative" samples and those of the second table as "differential" 
samples. 
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columns: 5W is the width of the innermost shell in 
/z"' Mpc; is the number of measures of the sample; 
(|cos 6\} is the mean of the |cos 6\; p is the characteris- 
tic parameter of Equation (j?]); and SNR is the signal to 
noise ratio computed with Equation (111. From simu- 



lations it has been found that results obtained with less 
than ~ 100 measures are not reliable, therefore, sam- 
ples with less than this number have been flagged with 
a question mark beside the value of the SNR. 

Errors in (|cos 6*1) are computed using the standard 
deviation of |cos 9\ resulting from 10000 Monte Carlo 
simulations with no signal (see below) corrected using 
Equations ( B1|[B3 1 assuming the following relation: 



cr(\cose\)(p + \\corr) _ (T(\co^e\)(p ^ 1; sim) 
o-(\cosB\)(p + \\thed) o-(|cose|>(f = \\theo) 



(12) 



where a-i^^o^Q\){p\thed) is computed using the the- 



oretical expressions described in Equations ( B1|B3 
0"(|cose|>(/' - 1; sim) is computed from the simulations 
and cr(|(.ose|>(p + 1; corf) is the final value used in Ta- 
bles gB 

Since p does not follow a Gaussian distribution, 
for this parameter we provide the confidence levels at 
Icr again correcting the theoretical values from Equa- 
tions ( |B5||B6| l with the confidence levels measured 
from the simulations with no signal. 

In Figures[3]and[4]are plotted the values of p (upper 
panels) and SNR (lower panels) as a function of the 
radius of the voids, for the cumulative and differen- 
tial samples, respectively. In the first figure the plotted 
radius is the minimum radius of each sample and in 
the second the central radius of each bin. In different 
colors are plotted the 10 shells widths that have been 
explored. 

The first result is that p < \ for most of the subsam- 
ples. This means that the direction of the spin of the 
galaxies tends to be parallel to the radial direction of 
the void. Despite of the large sample that we are using 
and the additional statistical correction, we find that 
the significance of the signal is not high most of the 
times and is dependent on the radius of the voids. The 
highest significance is reached when selecting galax- 
ies around voids, larger than 16/; ' Mpc in shells of 
3/i"' Mpc (\SNR\ > 3.6) using 179 galaxies, but the 
\SNR\ is higher than 3 increasing the width of the shell 
up to 7 Mpc and the sample size to 614 galaxies. 
Therefore, the further the galaxies are from the surface 
of the void, the lower is the signal although the in- 
crease in the size sample keeps the significance high. 



The relation between the strength of the alignment 
and the radius of the void is clearly shown in Fig- 
ure [4] where the differential samples are used. For 
voids smaller than ~ 15 Mpc these results are com- 
patible with a random distribution. It is for voids 
> 15 h^^ Mpc when it appears a signal of alignment, 
although given the smaller size of the samples the 
highest significance reached is 2.96 for voids with 
16/!"' Mpc < Rvoid ^ 17/!"' Mpc and a shell of 
6/i"' Mpc. 

Figure [6] shows the corrected histograms of 6 val- 
ues (see Equation |A5| l for the cases in which is reached 
the maximum S NR for the cumulative (left panel) and 
the differential samples (right panel). The continuous 
red line shows the analytical model described by Equa- 
tion (jT]) with the p values corresponding to these two 
maxima. 

To compute in a more robust way the significance 
level of our results we need to compare them with a 
control sample with no signal. To construct this con- 
trol sample, we have run 10000 Monte Carlo simula- 
tions in which the spin direction of the galaxies (de- 
termined by their position angles and axial ratios) has 
been shuffled so that each galaxy it is assigned the spin 
direction of any other galaxy randomly selected. This 
procedure has the advantage of ensuring the random- 
ness of the spin distribution and, therefore, the lack of 
any alignment signal, while using real data. 

For each simulation we have repeated the analysis 
performed in the real data using cumulative subsam- 
ples, and we have also computed 2 statistics used both 
in the real data and in the simulations. These statistics 
were computed as follows: for each bin in Rvoid, the 
mean (median) value of the SNR measured in the 10 
difl'erent shell widths was computed and the extreme 
value (i.e. with highest absolute value) of each simu- 
lation was kept. 

The significance level from each statistic is com- 
puted as 

1 - 2f(extreme{S NRs im} < extreme{5A^/?R£.fl/)), 

where f{extreme(SNRsim) < extieme{S NRueai]) is 
the fraction of simulations with extreme values of the 
mean (median) of SNR lower than the observed ones 
(-2.73 for the mean; -3.10 for the median). The fact 
of multiplying by two the observed fraction takes into 
account that we are considering only one side of the 
distribution of the statistics. 

The results from this analysis are summarized in 
Table [3] Using the mean we obtain a significance of 
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98.8% which improves slightly to 99.5% if the median 
is used instead. 

Although this test ensures the existence of a global 
signal, we have also checked that the increase of the 
SNR with the radius of the voids shown in Tables]?] 
and[5]is not just a consequence of the variation in the 
size of the samples. Figure ]5] shows the dependence of 
the SNR with the size of the samples (Ng). For clarity 
purposes we have restricted the analysis to the voids 
larger than 14 /i ' Mpc. Each line corresponds to dif- 
ferent limits in Rvoki and each point corresponds to a 
shell width. It can be seen that all subsamples show a 
similar trend with the width of the shell, showing the 
maximum value of SNR, in absolute value, at interme- 
diate shell widths. Also, for subsamples with similar 
sizes, the SNR shows variations depending on the size 
of the voids and shells, rejecting the hypothesis that 
the variations of the SNR are due only to variations of 
the size of the samples. 

Finally, to check the dependence of the signifi- 
cance of the signal with the distance of the galaxies 
to the surface of the void, we have constructed sub- 
samples containing galaxies closer than 5/?"' Mpc of 
the voids' surface and galaxies between 5 /i"' Mpc and 
10 /j"^ Mpc. With these two groups of galaxies we 
have repeated the analysis (only cumulative). In Ta- 
ble [6] are presented the results of this analysis which 
show that galaxies at distances larger than 5 /i ' Mpc 
do not present any significance aUgnment. 

6. Discussion 

We have analysed a large sample of galaxies around 
699 voids with radius larger than lO/i 'Mpc up to 
z - 0.12. We have found that for voids with ra- 



Criteria 



Ns 



miniiSNR)} < -2.73 58 
min{median(5A^/;)) < -3.10 23 



98. 



Table 3: Results from 10000 simulations with reshuf- 
fling of position angle and axial ratio between 
galaxies, showing the number of simulations pre- 
senting values for two statistics smaller than the 
observed values (extreme{(5A^7?))(Real Data)=-2.73; 
extreme{median(5A^7;))(Real data)=-3.I0). The last 
column shows the significance of the result, which 
takes into account that we are considering only one 
side of the distribution of the statistics. See text for 
more details. 



dius > 15 /z"' Mpc and within a shell not larger than 
~ 5/!"' Mpc, disk galaxies present a significant ten- 
dency to have their spin vectors aligned with the radial 
direction of the void. 

The maximum \SNR\ - 3.62 is measured for 
voids with Rvoid ^ 16 /j 'Mpc and a shell width 
of 3/i ' Mpc with a strength of the alignment p = 
0.664^Qy^i^. However, this value gives an overestima- 
tion of the real strength since has been selected as the 
best case out of many subsamples. 

In the next sections we compare our results with 
previous empirical works and with results from numer- 
ical simulations. 

6.1. Comparison with empirical works 



T06 and. 



From the observational point of view, 
more recently, ]S09| have performed a similar analysis 
to the one done here. T06 analysed 201 face-on and 
edge-on galaxies around voids with R > lO/i ' Mpc 
using data from the SDSS-DR3 and the 2dFRGS. They 
found a significant tendency of the spin of the galaxies 
to be in the direction perpendicular to the radial direc- 
tion of the void. More recently, |S09] using two samples 
of 578 and 258 galaxies from the SDSS-DR6 with sim- 
ilar selection criteria found no statistical evidence for 
departure from random orientations. 

Using the same criteria on the size of the voids 
(Rvoid > lO/j 'Mpc) and on the width of the shell 
(4 /i"' Mpc) as in those previous works, we find no sig- 
nificance alignment (p = 0.998; SNR = -0.15; see 
Table |4]l. The size of the sample used to establish this 
result is of 11060 galaxies, which after applying the 
correction factor of 0.6, means an equivalent size of 
6636 galaxies. This number is 8 times larger the size 



used by S09 



For a better comparison, we have computed the sig- 



nal of the alignment using criteria similar to that of S09 
regarding the selection of spiral galaxies (g - r < 0.6) 
and the definition of edge-on (b/a < 0.27) and face-on 
(b/a < 0.96) galaxie^ However, we keep our limit 
in Mr instead of using Mr > -21-1- Slogh as done 



by S09 because otherwise the final number of galaxies 
would be too small. After applying these criteria, we 
finished with a sample of 252 face-on and edge-on spi- 
ral galaxies. The value of p obtained with this sample 
is 0.993 with a SNR = 0.1 and, therefore, compatible 



'^We used an slightly different way to compute the axial ratio b/a 
compared with , 809, however, we do not consider this to have a sig- 
nificance effect on our comparison. 
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Fig. 3. — Upper panel: Values of the parameter p for subsamples of galaxies in shells of width SW and voids with 
radius larger than Rmui- Lower panel: Signal to noise of the alignment found for each subsample. This gives the 
significance of rejecting the null hyphotesis of not existence of any alignment. 
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0.569 
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0.520 



+ 0.003 
+ 0.004 
+ 0.005 
+ 0.006 
± 0.009 
+ 0.013 
+ 0.020 
+ 0.027 
± 0.039 



0.904! 
798+0 090 

^"-0.083 

924+0- 154 



-0.053 
-1.365 
-1.864 
-1.492 
-2.814 
-3.297 
-3.473 
-2.141 
-0.516 



30233 

19808 

12064 

6975 

3514 

1706 

723 

423 

189 



0.502 ± 
0.501 ± 
0.504 ± 
0.503 ± 
0.502 + 
0.504 ± 
0.495 ± 
0.491 ± 
0.502 ± 



0.002 
0.003 
0.003 
0.004 
0.006 
0.009 
0.014 
0.018 
0.027 



m 



-0.915 
-0.374 
-1.170 
-0.698 
-0.304 
-0.428 
0.350 
0.496 
-0.080 



Table 6: Strength of the alignment for galaxies within a shell up to 5/i ' Mpc (left) and for galaxies within a shell 
between 5 Mpc and 10 /i"' Mpc (right). Each line corresponds to samples of galaxies around voids with radius 
larger than Ry^"^^- Errors are computed using theoretical expressions in Appendixpl 



12 



2 - SW.1 Q Q SW.2 *— * SW.3 9-» SW=4 «-« SW=5 T-T SW.6 <— « SW.7 >-» SW=8 »-« SW=9 SW=1o| _ 




10 11 12 13 14 15 16 17 18 19 




-4- 



10 11 12 13 14 15 16 17 18 

Rvmd [h-^ Mpc] 



Fig. 4. — Upper panel: Values of the parameter p for subsamples of galaxies in shells of different widths SW (in 
Mpc) and voids in bins of 1 /i"' Mpc in radius. Lower panel: Signal to noise of the alignment found for each 
subsample. This gives the significance of rejecting the null hyphotesis of not existence of any alignment. 
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Fig. 5. — Variation of the SNR as a function of the number of galaxies A^^ for different samples. Each line corresponds 
to a selection in void's radius (differential in the upper panel and cumulative in the lower one) and each number 
indicates the width of the shell in /i ' Mpc. For practical purposes, only the samples of voids larger than 14 /i ' Mpc 
are shown. 
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Fig. 6. — Distribution of |cos 6\ after applying the statistical correction for the two subsamples that reach the highest 
\SNR\ in Table [4] (left) and Table [5] (right). The continuous line corresponds to the theoretical model described by 
Equation (|7]i for the measured value of p, shown in the upper left corners. 



with a random distribution. 

In relation to |T06| it is important to note that in 
that work, the signal presented corresponds to the peak 
of the signal found after exploring in shells of dilfer- 
ent widths. Consequently, to properly address the sig- 



nification of the signal in T06 we have reviewed the 
data used in that paper taking into account this ex- 
ploration. To do this, we have run simulations with 
similar number of galaxies around similar voids (R > 
10/;"' Mpc) searching for the maximum of the \SNR\ 
in shells of width from 3 /; ' Mpc to 7 /1 ' Mpc in steps 
of 0.1 /; ' Mpc. Then, it has been computed the frac- 
tion of simulations with max(\S NR\) larger than the 
one found in T06 when only SDSS data was used 
(max(S NRrob-^sDss) - After doing our analy- 

sis we found that ~ 16% of the simulations showed 
max(\S NR\) > 2, decreasing the significance of the re- 
sult oflTOe'to ~ 84%. 

Another observational work on alignment of galax- 
ies with respect to the local large scale structure is 
that by Lee & Erdogdul ( 2007 1. In this paper, the au- 
thors computed the shear tidal tensor in the position of 
each galaxy and measured the angle between each of 
the principal axes and the spin direction of the galaxy. 
They found a significant alignment (> 6cr) between 
the direction of the spin and the intermediate princi- 
pal axis of the shear tidal tensorp] Nevertheless, the 



'^In the process of conducting this analysis we noted some duplica- 
tions in a few galaxies that make the measured SNR decreases from 
2.4 to 2. 

'''The authors obtained a value of c = 0.084±().014 which con'esponds 



comparison with our results is not direct since we do 
not use a direct measure of the orientation of the shear 
tidal tensor in the position of each galaxy and the ra- 
dial direction of the voids can be considered only as a 
statistical proxy for the direction of the major princi- 
pal axis of the shear tidal tensor. Given the dilferences 
in methodology, a meaningful comparison of the re- 
sults from both works would need an analysis that it is 
beyond the scope of this paper. 

6.2. Comparison with numerical simulations 

Another way to study the alignment of galaxies 
with their local large scale environment is through nu- 
merical simulations ( |Porciani et al.| [2002a|b[ [Navarro 
et al. 2004; Baili n & S teinmetz' 2005 ; Bailin et al. 
2005, ,Altay et al. 2006 



et al 



2008 



2010 



2007 



Aragdn-Ca 



Patiri et al.,. 2006b JBrunino 



vo et al. 



2007^ 



Cuesta et al. 



Zhang et al. 2009; H ahn et al.|20 



Bettetal. 



Since the behaviour of the halos can be de- 
pendent on the environment or the large scale struc- 
ture in which they reside, to perform a meaningful 
comparison we have focused on the analysis done by 
Patiri et al.| ( |2006b| l, |Brumno et al.| ( [2007l ) and |Cuesta| 
et al. ( 2008| l in which it was studied the orientation 
of dark matter halos around cosmic voids using dif- 
ferent cosmological simulations. The criteria imposed 
to the dark matter halos and the procedure to detect 
voids tried to match the criteria used in iT06i All three 
works found that the minor axis and the major axis 



top= 1.057 ±0.012 
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of the halos have significant tendencies to be ahgned 
with the radial and the perpendicular directions, re- 
spectively. The results regarding the orientation of 
the angular momentum of the halos were less clear 



Patiri et al. ( 2006b i did not find any particular orienta- 



tion for the angular momentum of the halos. |Brunino| 
|et al.| ( |2007| l also did not found any particular align- 
ment for the angular momentum in their full sample of 
halos although they detected a tendency for those halos 
with a disc-dominated galaxy to have their angular mo- 
mentum perpendicular to the radial direction. Finally, 
|Cuesta et al.| ([2008 ) measured a significant (> Icr) ten- 
dency of the spin of the dark matter halos to lie in 
the plane perpendicular to the radial direction. How- 
ever, the same authors found that the strength of the 
alignment is mainly produced by the outer regions of 
the DM halos and this would explain the discrepancies 
with |Brunino et al.] ( |2007| l were the inner regions of the 
DM were used to measure the alignment. 

How these results relate with ours is not straightfor- 
ward since we observe the collapsed baryonic matter 
and they studied the dark matter or the non-collapsed 
baryonic matter We can only point out the fact that 
the alignment that we find in the galaxies is shared by 
the minor axis of the halos studied in the simulations, 
either dark matter or gas halos. This is suggestive to 
an interaction between the galaxy and the hosting ha- 
los around it (either of dark or baryonic matter) lead- 
ing to a tendency of the minor axis of the galaxy (and, 
therefore, its angular momentum) to be aligned with 
the minor axis of the halo's matter distribution. 

7. Summary 

Analysing a volume of ~ 27 x 10^ (/i"' Mpc)-' from 
the SDSS-DR7 we have searched for cosmic voids de- 
void of galaxies brighter than M,- - 5 log/i - -20.17 
and with Rvoid > 10/; ' Mpc. We have found 699 non 
overlapping voids for which we provide positions and 
sizes. 

We have used this catalog of voids to search for disk 
galaxies around them and study the alignment between 
the direction of the angular momentum of these galax- 
ies and the radial direction with respect to the center of 
the voids. 

We have included two improvements with respect 
to previous similar works. 

First, we have used an updated version of the SDSS 
spectroscopic catalog (data release 7) and we have 
combined it with the visual morphological classifica- 



tion from the Galaxy Zoo project to get a reliable sam- 
ple of disk galaxies. 

Second and more important, we have introduced a 
statistical procedure that has allowed us to overcome 
the problem of the indetermination of the real incli- 
nation of galaxies computed from their apparent axial 
ratio. We have performed extensive Monte Carlo sim- 
ulations to check the validity of this procedure. We 
show that the procedure recover the real signal without 
practically any bias and its power in terms of capacity 
to reject the null hypothesis it is equivalent to the case 
in which it is used a sample with complete knowledge 
of the real direction of the spin of the galaxies using 
60% the amount of galaxies. In comparison with the 
common procedure of selecting only edge-on and face- 
on galaxies, this procedure means an increase of about 
a factor 3 in the amount of measurements used in the 
analysis of the alignment. 

These improvements have allowed us to detect a 
statistically significant (> 98.8%) tendency of galax- 
ies around very large voids {R > 15/i"' Mpc) to have 
their angular momentum align with the radial direction 
of the voids. However, for smaller voids this tendency 
disappears and the results are consistent with no spe- 
cial alignment. 

We have also found that the strength of the align- 
ment depends on the distance of the galaxies to the 
surface of the voids and for galaxies further than 
~ 5 Mpc the distribution of the alignments is com- 
patible with a random distribution independent of the 
size of the voids. 

Previous similar works found opposite ahgnment 
( |T06 I) or no alignment (S09 ). However, these works 
used too few galaxies around voids with R > 10 /i ' Mpc 
which, according to our work, could mask the signal. 
In fact, using the same criteria for the size of the voids 
and the width of the shells as in those works, our data 
is compatible with a random distribution of spins with- 
out any particular alignment, as found by 809' 

The comparison with the results from cosmolog- 
ical simulations points to a possible connection be- 
tween the alignment of the halos (of dark matter and 
non-collapse baryonic matter) and that of the galax- 
ies which could explain the similar orientation of both 
components observed in the simulations and in our 
work, respectively. 
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A. Statistical computation of P(|cos 0\) with full information 

In this Appendix we describe in detail the statistical procedure that has been used to compute the corrected value of 
p using all the disk galaxies, independently of their inclination. We also show its validity and robustness using Monte 
Carlo simulations. 



A. 1 . Mathematical justification of the procedure 

In the analysis of the alignment of galaxies, one of the main sources of uncertainty is the indetermination in the 
inclination of the plane of a disk galaxy with respect to the line of sight using exclusively the observed axial ratio 
b/a of the galaxies (see Equation[T]i. In other words, if a galaxy is divided in two halves separated by its major axis, it 
is not possible to know which of the two halves is the closest to the observer]^ This indetermination is negligible for 
edge-on (b/a ~ 0) and face-on galaxies (b/a ~ 1) but using only these galaxies reduces the sample size to ~ 1/5 of 
the original one. As we increase the range of allowed values of b/a, the increasing uncertainties in ^ will result in an 
increasing degradation of any existing alignment but the statistics improve. The question is whether this improvement 
of the statistics can compensate for the increasing degradation of any possible signal. The answer is "yes". Choosing 
always the plus sign in the computation of ^ , or the minus sign, or any random assignment of signs, leads to the same 
statistical results (i.e. they are equally powerful tests), which are better than those obtained with any limitation of the 
range of possible values of b/a. However, the estimate of the alignment obtained in this manner is biased towards 
smaller values (the strength of the alignment is given by (I - p) ^ -3c/4 for weak alignments). This would not be 
much of a problem, because one may calibrate the procedure using Monte Carlo simulations and then correct for the 
biasing. In this manner we have found (see Table]?]): 

1 - p^^ 0.6(1 -P) (Al) 



where p corresponds to the real alignment and p+ is the value obtained using the plus sign for f . 

The main problem with the use of p+, or any other sign assignment, is that it introduces an artificial randomness that 
increases the scatter of the estimates. Using p+, we assign the correct sign to half of the galaxies, on average, while 
the other half gets the wrong sign, but the exact number of galaxies getting the correct sign fluctuates from sample to 
sample (with variance Ng/4, being A^g the size of the sample) resulting in an increased error Furthermore, since we do 
not take into account the other possible sign assignment, we do not know how large is the degradation of the alignment 
implied by those galaxies that get the wrong sign. 

To avoid these problems we propose a method that uses all the information in the data and does not introduce 
artificial randomness. To this end, we consider the two possible values of associated with every galaxy (one value 
for each possible sign of and assume that only half of the values of 6 falling in a given range are correct while the 
other half is incorrect. The correct values for the latter half of galaxies would be the conjugate of 6, 0' , corresponding 
to the value of 6 using the opposite sign for Thus, if the actual probabihty distribution for 6 werej^ 

P(cos e, p) = r — , , (A2) 

' (1 +(/?2 - l)cOs2 6')3/2) ^ ' 

the probability distribution that would be inferred from the 2Ng values of 9 treating them as if they were independent, 
f (cos ff), would be given, for the j - th bin, by: 



1- 1 , 

P(cos Oj) = -P(cos 0j, P)+Yi}^ Z')' (A3) 



/■=i 



where / is the number of Oj values in the j - th bin. 



The use of kinematic information or tlie presence of dust lanes can lielp to break this indetermination. however, in most of the cases this is 
information is not accessible. 

'We are only interested in the direction of the alignment and therefore the analysis can be restricted ioO <6 < nl2 and cos 6 = |cos 6\. 
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This formula expresses the fact that with probability 1/2, the probability density in the bin centered in Oj is given 
by the real distribution P evaluated at 6j (correct sign assignment), while with probability 1/2, the probability density 
at 6j is the averaged of the value of P over the conjugate values {0'j{i)) of the I values of 9 falling in bin j. 

So, the factor: 

P(cos Oj) 

iP(cos Oj, p)+^, Picos 0'jii), p) 



2(cos Oj) ^ \ '\ ^ (A4) 



is an estimate of the ratio between the actual distribution, P(cos 6j), and the first estimate, P(cos Oj). 
Therefore, we have for the estimate of P (that we denote by Pc): 

P,.(cos dj) = P(cos 0j)Q{cos Oj) (A5) 



When the alignment is very strong, the assumption that the two values 6, 0' of a conjugate couple have the same 
probability can no longer be mantained. Instead, we should used: 

P(cos 0) 

Prob(6l) = (A6) 

P(cos 6) + P(cos 6') 

Pmh(e') = ^ (A7) 

^ ' P(cos6) + P{cos6') ^ ' 

and modify the definition of Q consequently. However, this complication of the method is not worthy to our 
purpose. In fact, from Table]?] we can see that even for considerable alignment strengths, the bias implied by neglecting 
this last refinement is small, and can be corrected by the following expression: 

PAb = 1 - (1 + 0.1(1 - p,f ){\ - Pc) (A8) 



where pc is the value obtained using Equation ( A5 i, and pdt is the debiased value. 

From Table |7] we can also see that the relative error, (Tp/\l - p\, is always larger for p+ or than for pc- For 
weak alignments (|1 - p\ < 0.1) the former is ~ 20% larger than the latter, while for larger alignments the difference 
diminishes. 

Finally, it must be noticed that the method that we have just described does not depend on the form of P(cos ff). 

A.2. Description of the procedure 

In this appendix we describe the actual implementation of the method presented above. 
The procedure is as follows: 



1 . For each galaxy, we compute the two possible values of cos corresponding to the two alternatives signs of ^ 
and hence the two possible spin orientations. 

2. Then, we construct a normalized histogram assuming the two values of cos of each galaxy as independent 
values. The normalization is done dividing each bin by 2 times the total number of galaxies of the sample (Ng) 
and by the width of the bins. We call this non-corrected histogram f (cos 0). 

3. Next, in each bin centered in cos0y, we compute the value of the corrected histogram Pc{cos0) using Equa- 
tions ( |A5] l and (|A4|: 



Pcicos 0j) = P(cos 0j) 2(cos 0j) (A9) 

remembering that 
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p(cosej,p) 

ip(cos 6j, p)+^, Picas e'j(i), p) 



and that 0'j(i) are the conjugate values of for those galaxies with 6j(i) within the interval | 6j(i) - Oj \< A6j/2, 1 
is the total number of values within the bin, and 

P(cose,p)^ ^ P (AlO) 

The final corrected value of p is computed numerically using its relationship with (cos 6} from Equation ([8]l, which 
given the distribution P^. (cos Oj) as a discrete distribution can be expressed as: 

Z"j=lPc(cOS0j) l+p 

with n the total number of bins in which the distribution is divided. 

A.3. Robustness of the statistical correction 

To check the robustness of the statistical correction we have performed a series of Monte Carlo simulations. In 
these simulations, we use samples of fake galaxies in the position of the real ones but with spin directions assigned 



randomly following a p - distribution (Equations A2 1 with a given pinput- Then, the samples of fake galaxies are 
analysed in the same manner of the real galaxies and a final Poutput value is obtained. 

We have run 2 sets of simulations using two samples with different number of galaxies {Ng) to check the robustness 
of the procedure also as a function of the sample size. These samples correspond to galaxies in shells of 4/i ' Mpc 
and Rvoid > 10 /i"' Mpc (Sample A, following the usual criteria used in previous works) and to galaxies in shells 
of 3 /i"' Mpc and RvoU > 16 h^^ Mpc (Sample B, corresponding to our maximum SNR). For each sample we have 
run 1000 Monte Carlo realizations with 7 diff'erent initial distributions of (|cos0|) described by their corresponding p 
values (pinput)- These values covered the typical values of p that we have found in our analysis. 

Table |7] shows the results of this analysis. For each subset of 1000 realizations we give the size of the sample 
A^, the input value Pi„p„i, the mean value of p obtained if a fixed sign for ^ is used (p+ and p_, for plus and minor 
sign, respectively) and the mean value of p when applying our stastical correction, p output- The uncertainties shown 
correspond to Icr of the distribution of the single values in the 1000 realizations. 

We found that for most of the cases the statistically corrected value is within Icr of the input value showing the high 
accuracy of the procedure, especially when comparing with the cases in which a fixed sign is used. 

The results of these simulations have been used to compute the "effective size" of the initial sample. This effective 
size is defined as the size that a sample with complete knowledge of the real signs of ^ for each galaxy should have to 
show the same uncertainties that we find in our simulations. On what follows, it is described how we have computed 
the correction factor to be applied to our samples to obtain their effective sizes. 

It can be proved theoretically that the value of the standard deviation of (|cos G\), cr(\co^e\), for the case in which there 
is no preferential alignment (p - 1), is: 

o"(|cosfl|> = ; , (A12) 

VT2A^ 

where is the total number of galaxies used to compute (|cos0|). However, this theoretical expression assumes 
the full knowledge of the values of 6 for all the galaxies, while empirically we do not have such full information 
because of the indetermination in the sign of f . Therefore, we have compared the standard deviation obtained from the 
simulations with different values of A^, with the theoretical value. From this comparison, we have obtained a correction 
factor to be applied to the total number of galaxies equal to 0.6. This means that our statistical approximation carried 
an uncertainty that is equivalent to the uncertainty of having ~ 60% of the galaxies with full information. 
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Sample pi„pui p+ P- {-'output 













iP-} 


o-p_ 


iP out pill} 


Foul put 


A 


11060 


0.50 


0.664 


0.008 


0.665 


0.008 


0.437 


0.010 


A 


11060 


0.75 


0.841 


0.010 


0.842 


0.009 


0.746 


0.012 


A 


11060 


0.90 


0.938 


0.010 


0.938 


0.010 


0.900 


0.013 


A 


11060 


1.00 


1.000 


0.011 


1.000 


0.011 


1.000 


0.014 


A 


11060 


1.10 


1.060 


0.012 


1.060 


0.012 


1.102 


0.016 


A 


11060 


1.25 


1.143 


0.012 


1.142 


0.012 


1.256 


0.018 


A 


11060 


1.50 


1.272 


0.013 


1.272 


0.014 


1.542 


0.025 


B 


179 


0.50 


0.642 


0.060 


0.644 


0.059 


0.486 


0.059 


B 


179 


0.75 


0.833 


0.071 


0.830 


0.075 


0.753 


0.082 


B 


179 


0.90 


0.938 


0.081 


0.938 


0.085 


0.907 


0.101 


B 


179 


1.00 


0.999 


0.088 


1.004 


0.086 


1.004 


0.107 


B 


179 


1.10 


1.073 


0.091 


1.072 


0.091 


1.115 


0.121 


B 


179 


1.25 


1.162 


0.101 


1.155 


0.100 


1.258 


0.140 


B 


179 


1.50 


1.305 


0.113 


1.313 


0.110 


1.573 


0.561 



Table 7: Results of several simulations to test the validity and robustness of our statistical correction. Two samples 
with different number of galaxies are shown: Sample A is made of galaxies in shells of Ah^^ Mpc around voids with 
R > 10 /i ' Mpc and Sample B is made of galaxies in shells of 3 Mpc around voids with R > 16 /z ' Mpc. Each row 
corresponds to a set of 1000 realizations in which to each real galaxy a synthetic spin vector was assigned following 
the theoretical distribution given by Equation (|7ji with a p - pmpui- p+ and p^ are the values of p obtained when fixing 
the sign of ^. Poutput is the final value after applying the statistical correction. For each parameter Pouipm)^ the 

mean and the standard deviation of the 1000 realizations are shown. See text for more details. 



B. Computation of the uncertainties in <|cos 6\} and p 

In this section we present the expressions used to compute the uncertainties in {|cos0|) and p, in the general case. 

The standard deviation of (|cos 6\}, cr^ij-osei), is computed as a-\cose\l -^JN^^ where Ng is the total number of galaxies. 
o"|cose| is the root mean square of |cos0| for the distribution given by Equation [7] Computing o-\cosB\ analytically, we 
find the following expressions for cr^^^g^g^') depending on the value of p - (|cos - 1, : 




ln(p + Jp^ - 1) - - T-^; P>1 

^ - 1 1 + 



(Bl) 
(B2) 



arcsin(Jl-p2)--^; p<l (B3) 
' 1 + / 



When using the method described in Appendix A. 3 we have some uncertainty in the direction of the spin, but we 
find that the errors are well described by the above expressions using 0.6 times the number of galaxies in the place of 
Ng (see Table|7]). 

Since the distribution of p is not Gaussian, we can compute the value p and the Umits of the Icr confidence interval 
(p-cT, P+it) with the next expressions: 

p = ^- 1 (B4) 

^ <|COS0|) 

= n ^ 1 

<|cos e\} + crQco!,0\) 
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